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Abstract 

At present the color-octet mechanism is still an important and debatable part in the non- 
relativistic QCD(NRQCD). We find in this work that the polarized double charmonium production 
at the LHC may pose a stringent test on the charmonium production mechanism. Result shows 
that the transverse momentum(pr) scaling behaviors of double J/ip differential cross sections in 
color-singlet and -octet production mechanisms deviate distinctively from each other while pT is 
larger than 7 GeV. In color-octet mechanism, the two J ftps in one pair are mostly transversely 
polarized whenpT ^ 2mc, as expected from the fragmentation limit point of view. In color-singlet 
mechanism, there is about one half of the charmonium pairs with at least one J/^ being 
longitudinally polarized at moderate transverse momentum. The energy dependence of the 
polarized J/ip pair production is found to be weak, and this process is found to be experimentally 
attainable in the early phase of the LHC operation. 



PACS numbers: 12.38.Bx, 13.85.Fb, 14.40. Lb. 



I. INTRODUCTION 



nrqcd[i| has now become a basic theory to parameterize non-perturbative contributions 
in heavy quarkonium production and decays by color-singlet (CS) and color-octet (CO) matrix 
elements. It defines the "velocity scaling rules" to order various matrix elements by the 
relative velocity of heavy quarks in the rest frame of heavy quarkonium. Although NRQCD 
has already achieved a lot in the study of heavy quarkonium physics, there are still some 
unsettled issues remaining, especially, in aspect of production mechanism {2!, for instance, 
the role of color-octet in charmonium production. 

The color-octet scheme is a central essence of NRQCD, which supplies a way to partly 
understand the prompt J/ip and ip' surplus production at the Fermilab Tevatron [3], and 
supplies a consistent differential cross section spectrum, versus the transverse momentum 
Pti with experimental results j4|. However, there are still something unclear in this scheme, 
which ask for further effort. In electron-positron collision, the former crises in J /if) exclusive 
js-?! and inclusive production js], [ ol are somehow going to an end with the advent of the 
next-to- leading order (NLO) QCD [10|, lHj and relativistic correction results 12|. In char- 
monium photoproduction, the recent NLO QCD calculation tells that to reproduce the full 
HERA experimental results, the color-singlet contribution alone seems not enough jl3 |. In 
charmonium hadropro duct ion, recent computations of NLO QCD corrections [14] signifi- 
cantly enhance the J/ip color-singlet yield in large transverse momentum region and alter 
its polarization prediction in LO calculation. Hence, the NLO J /ip hadroproduction result 
in CS scheme greatly minimize the contribution through color-octet mechanism, which in 
some sense agrees with the hypotheses given in Ref. |15|. For more detailed discussion on 
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these issues, readers can refer to recent review articles, e.g. 

To further investigate charmonium production mechanism keeps on being an urgent and 
important task in the study of quarkonium physics. The running of the LHC supplies a 
great opportunity to this aim. With a luminosity of about 10'^^ ~ 10'^^cm~^s~^ and center 
of mass energy 7 ~ 14 TeV, LHC will produce copious data of charmonium inclusive and 
exclusive production, which can in principle answer the question: to what extent the color- 
octet mechanism plays a role in tp production? As noted by Cho and Wise that the direct 
ip production via color-octet mechanism should be mostly transversely polarized at leading 



order in a<(, which may stand as a unique test to the quarkonium production mechanism 



19(1 . In this work, we consider the process of polarized J/ip pair production at LHC, and 



suggest to test the charmonium production mechanism through it. The same process but in 
unpolarized case was proposed and calculated by Barger et al. 20 1 and Qiao once for 



Tevatron experiment. 

The paper is organized as follows: In section II, we set up the formalism for the compu- 
tation of polarized if) hadropro duct ion, in both color-singlet and -octet schemes. In section 
III, the numerical results for the LHC first-year and later runs are provided. The section 
IV is remained for a brief summary and conclusions. For the convenience of readers, the 
relevant analytic expressions are given in the Appendix. 



II. FORMALISM 

A. Color-Singlet Scheme 

In color-singlet model, the partonic subprocesses start at order af, which include g + g ^ 
J/ip + J/ip and q + q ^ J/ip + J/ip. Intuitively, the latter, the quark-antiquark annihilation 
process, is possibly negligible in comparison with the former since the LHC is a proton-proton 
collision machine. To confirm this point, the quark-antiquark process is also evaluated in 
this work and the numerical result is given in Table HI Hence, in the following analysis we 
will mainly focus on the gluon-gluon process as shown in Figure [H and the analytic results 




FIG. 1: Typical Feynman diagrams of J ftp pair production in pp collision at leading order in 
color-singlet scheme. 



of it are presented in the Appendix. 

For J/ip{^Si) color-singlet hadronization, we employ the following commonly used pro- 
jection operator 
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v{p,)u{Pc) ^ fji^ (f + 2m,) X ( ^^^/^(O) ) ® ( ^ ) , (1) 



where e^^^ is the J/?/' polarization vector with P ■ e = 0, 1, stands for the unit color matrix, 
and A*", = 3. The nonperturbative parameter, '?/'j/^(0), is the Schrodinger wave function of c 
and c system at the origin for J/V^, and |V^j/^(0)| = ^\R{{S)\^ j ^-n with P(0) being the radial 
wave function at the origin. In our calculation, the non-relativistic relation mjj^ = 2mc is 
adopted, and hereafter, for simplicity mj/^ will be denoted by m. 

The differential cross section for J/ip pair hadropro duct ion reads as 

da , _ , /",,„. ^ ^_ da 



dpT 



dyidy2fa/p{xa)fb/p{xb)'2pTXaXb-^{a + SJ/z/;) , (2) 



where fa/p and fb/p denote the parton densities; yi, y2 are the rapidity of the two produced 
J/ips. The partonic scattering process, the gluon-gluon to polarized J/tp pair can be 
calculated in the standard method. To manipulate the trace and matrix-element-square for 
those diagrams shown in Figure [H the computer algebra system MATHEMATICA is em- 
ployed with the help of the package FEYNCALC and FEYNARTS [22]. The lengthy expres- 
sions for different J/ ip pair polarizations are presented in the Appendix, for the convenience 
of independent comparison and physical simulation use. Sum ming over all polarizations, we 
find an agreement with the unpolarized cross section given in 21 1. 

In experiment, the J/ip polarizations can be determined via measuring the normalized 
angular distribution I{cos6)of ^~ pair in J/ip decays. The correlation of distribution 



with polarization reads 23|, |24| 



/(cos^) = — — + a cos^^). (3) 
l[a + 6) 

Here, 6 is the angle between the /i"*" direction in the J /ip rest frame and the J /ip direction in 
pp center-of-mass frame; a is a convenient measure of the polarization, which is defined as 
a = (cTy — 2aiP)/{aT + ^ai) where ar and a^ are transverse and longitudinal components of 
the cross section, respectively. For detailed derivation of the relation between polarization 



measure a and decay distribution, see, e.g., the Appendix A of Ref. 



24j . From the definition. 




FIG. 2: Typical Feynman diagrams of J/^ pair production in pp collision at leading order in 
color-octet scheme. 

obviously, a — 1 corresponds to the transverse polarization, and a = — 1 to the longitudinal 
polarization. 

B. Color-Octet Scheme 

In NRQCD, the short-distant contribution can be expanded by v. Dynamical gluons enter 
into Fock state, and combine with heavy quark pair to form Color octet states. Therefore, 
wave functions of J/ip is represented by 

I JM = 0{l)\c-c['s['^]) + 0{v)\c-c['PP]g) + 0{v')\c-c['s[''%g) + 0{v')\c-c[' 3^%) + ... (4) 

In the CO scheme, at the leading order of v, three different CO state components 
\cc[^Pj''']g), \cc[^S[^'^^]gg) and \cc[^ S^'^'lg) may contribute to the double J/i/j production. 
In practice, for the sake of simplicity for further analysis, we classify the intermediate 
Fock states into two clusters: one contains \cc[^S[^'^^]gg), the other includes \cc[^SQ^^]g) 
and Ice fP!j%). In this sense, the double J/ip Fock states may have following different 
combinations: 

" V ' 

Parti 

+ \cc['sP]){\cc['PP]g) + \cc['si%)) + \cc['S[%g)\cc['S[%g) 

^ V ' ^ V ' 

Part2 Parts 

+ \ccM%g){\cc['PP]g) + \cc['S^%)) 
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+ i\cc?PP]9) + \ctS^%)){\cc['PP]g) + \cc['S^%)) (5) 



Parts 



In (E]) the leading one indicates the pure color-singlet state, which is discussed in the pre- 
ceding subsection. The CO states involved components are then divided into five sectors. In 
small transverse momentum pt region, where the NRQCD factorization theorem does not 
work well, the five CO state related processes would contribute less to the double charmo- 
nium production than from the pure CS process due to a suppression in matrix elements 
of or v^. In large pt region, however, the CO contribution, especially from |ccpS'|^^](y'(y') 
may exceed what from the CS, because the relative smallness of CO matrix elements may 
be compensated by the enhancement of the large propagators in corresponding processes, 
which is similar to the case of inclusive ip production in large transverse momentum region at 



the Fermilab Tevatron jj, l25|]. In infinite transverse momentum limit, it is well-known that 
the process of gluon fragmenting into \cc[^s[^^]gg) will dominate over others in hadropro- 
duction. Hence, one may still expect here that in the double J/ip exclusive hadroproduction 
in the large transverse momentum region, the bi-gluon fragmentation would be the most 
important production mechanism. However, in moderate transverse momentum case as in 
LHC experiment, to minimize the uncertainties induced by the fragmentation mechanism 
we calculate the complete Feynman diagrams for CO processes instead of fragmentation 
simplification as shown in Figure O 

For processes involving the CO components, that is parti to part^ in (|5]), those contribute 
double J/tp production via intermediate states \cc[^S^^]g) and \cc[^Pf'^]g) are obviously 
negligible in comparison with the CS process as shown in Figure [H Moreover, the part^ 
involved processes may contribute less than what from parti processes by f ^, an order. For 
Parti contribution, our explicit numerical calculation shows that it is smaller than what 
from the CS by about two orders of magnitude in all px region. In all, the meaningful CO 
contribution mainly comes from the part^, involved processes, and which are what we are 
going to evaluate. 

For the J/ip intermediate Fock state \cc^Sf'\gg), the projection operator is 

^8(0) 



v{p,) u{Vc) -j==T^cc fji^ ( r + 2m,) , (6) 



\Rm?=''-<0'l\^Si)> . (7) 



Except for the difference in CO and CS non-perturbative matrix element projections, tlie 
perturbative calculations of Feynman diagrams for both CS and CO are similar. 



III. NUMERICAL RESULTS 



In numerical calculation, we enforce the LHC experimental condition, the pseudorapidity 
cut \ri{J/il))\ < 2.2, on the produced charmonium pair. The cross sections and transverse 
momentum spectra are evaluated under three conditions: the first one with center-of-mass 
energy \fS = 7 TeV, the second one with center-of-mass energy y/S = 10 TeV, the third one 



with center-of-mass energy \/S = 14 TeV, which correspond to t 



le experimental situations 



of initial and following runs for hadron study at the LHC [26|]. The planned collisions 
happen at the LHC with luminosity of ~ 10^^ cm~^s~^ in first two cases, and ~ 10^'^cm~^s~^ 
in the third situation. In one year running, lO^s effective time, the corresponding integrated 
luminosities are about and 10(/6^^), respectively. The input parameters take the 



values 
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rric = 1.5 GeV, |R(0)|' = 0.8 GeV^ < C8^'''(^Si) >= 0.012 GeV . (8) 

The typical energy scale in partonic interaction is set to be at mx = \/ + p^, and hence 
the strong coupling is running with transverse momentum. 

With the above formulas and inputs, one can readily obtain the polarized J/ip pair 
production rate at the LHC. In our numerical calculation, the parton distribution func- 



tions(PDFs) of CTEQ5L 271] is used, and both renormalization scale in strong coupling as 
and factorization scale in PDFs are set to be m^. The numerical results for integrated cross 
section a{p p — t- J/ipJ/ip) with different px lower bounds are presented in Table I for the 
CS production scheme, and Table II for the CO scheme, where the the branching fraction 
of B{J/ip — )■ = 0.0597 is taken into account. 

The spectra of double- J/?/' exclusive production as function of transverse momentum 
Pt are illustrated in Figures [3] to From these figures we see that the conventional CS 
production scheme dominates over the CO one in relatively low-px region, some pt < 7 
GeV. 



TABLE I: The integrated cross sections of J /if) pair production in color-singlet model under various 
transverse momentum lower cuts at the center-of-mass energy ^/S = 7 TeV, 10 TeV and 14 TeV. 
Here, _L± represents for the situation of both J/'i/^s being transversely polarized, |{|{ for both J/tps 
being longitudinally polarized, ||_L for one J/-^ being longitudinally polarized and the other being 
transversely polarized. The totgg and totqq denote the polarization-summed cross sections of gluon- 
gluon fusion and quark-antiquark annihilation processes, respectively. 
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TABLE IL The integrated cross sections of J/ij; pair production in color-octet model with various 
transverse momentum lower cuts at the center-of-mass energy yfS = 7 TeV, 10 TeV and 14 TeV. 
Here, J— Lgs represents for the situation of both J/^jjs being transversely polarized, for both 
J/V'S being longitudinally polarized, H-Lgs for one J/V' being longitudinally polarized and the other 
being transversely polarized, and the iotss means the polarization-summed CO cross section in 
gluon-gluon fusion. The totis in last row represents for the double J/ip yields from CS + CO 
production scheme for reference. 
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IV. SUMMARY AND CONCLUSIONS 

We have in this work evaluated the polarized J/%1) pair hadroproduction rate, directly 
and exclusively, at the LHC in color-singlet and -octet schemes. In color-octet scheme, Fock 
state \cc}iS'f^]gg)\ccf'S'f^]gg) involved process dominates over others, and it yields mainly 
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FIG. 3: The differential cross-section of J/ip pair production versus pT at the LHC in color-singlet 
scheme. Lines from top to bottom, i.e. a ,b and c, represent yields from color-singlet in _L_L, ||||, 
and ||_L cases, respectively. The diagrams from left to right represent for cases o^^/S = 7 TeV, 
= 10 TeV, and VS = U TeV, respectively. 

the transversely polarized J/ips. Other CO related processes, either single or double CO 
processes, contribute significantly less than those from CS and parti involved processes in 
whole transverse momentum region, therefore can be neglected. Result shows that in low 
transverse momentum region the double J/ip yields are governed by the CS process, while 
the CO process may dominate as pt > 7 GeV. The indirect J/ ip production from the Xcj 
and ip' feed down, according to the analysis of {20], may contribute about a factor of 2 to the 
transverse distribution spectrum of prompt J/ip production, while these inclusive yields can 
be removed in experiment in principle and hence are not taken into account in our analysis. 

The calculation performed in this work is at the leading order of strong coupling ag. Due 
to the relatively high interaction scale in this process, higher order QCD corrections should 
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FIG. 4: The differential cross-section of J/ip pair production versus px at the LHC in color-octet 
scheme. Lines from top to bottom, i.e. a ,b and c, represent yields from color-octet in _L_L, ||_L, 
and nil cases, respectively. The diagrams from left to right represent for cases of s/S = 7 TeV, 
^/S = 10 TeV, and ^/S = 14 TeV, respectively. 

be small. The " smallness" here means that it may not be as big as NLO corrections to some 
inclusive charmonium production processes recently found. Nevertheless, this is just a naive 
speculation. For exclusive process, higher order corrections, the virtual corrections, do not 
change the relative polarization ratio, while is different from NLO correction to the inclusive 
process where the hard real correction may alter the final state polarization significantly. 
Though higher order relativistic correction is generally large in charmoniun system, while 
for J/tp pair production it will be doubly suppressed. It is worthy of mentioning that in 
order to keep the exclusivity of double J/ip production and reduce the final state interaction 
effects, to enforce the transverse momentum cut veto would be a necessary choice, however 
which may also suppress some useful signals. 
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FIG. 5: The differential cross-section of J/ip pair production versus pr at the LHC. Lines a and b, 
represent yields from color-octet and color-singlet in unpolarized case, respectively. The diagrams 
from left to right represent for cases o^^/s = 7 TeV, VS=10 TeV, and = 14 TeV, respectively. 

The integrated cross sections of double J/ip production with different polarizations are 
found to be measurable theoretically with a lower transverse momentum cut of 7 GeV 
in the first-year run of LHC at colliding energy of 7 TeV and luminosity of Ifb^^ or so, 
where about one hundred of di-muons from J/ip decays will be produced and the total 
yields from CS and CO are comparable. For lower transverse momentum cut higher than 
7 GeV, the double J/ip yield may mainly come from the color-octet scheme, under the 
circumstance of employing the nowadays prevailing magnitude of color-octet matrix element. 
At colliding energy of 7 TeV, our calculation result shows that the data with at least one 
J/tjj in a. pair being longitudinally polarized are nearly fifty percent of the total yield with 
transverse momentum integration lower cut of 5 GeV of produced J/ip. Different from the 
inclusive production, where experimental data show that J/ips tend to be unpolarized or 



even longitudinally polarized with the increase of transverse momentum, in the exclusive pair 
production both color-singlet and -octet schemes exhibit that the produced J/tp data should 
be always dominated by the transversely polarized yields. The above conclusions might be 
testified in future LHC experiment, and we expect the experiment measurement on exclusive 
double J /ip production may tell us more on the charmonium production mechanism. 
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Appendix 

The polarized differential cross sections in color-singlet model. Here, the symbol |||| 
represents for the situation of both J/ip being longitudinally polarized; ||_L for one J/iIj 
being longitudinally polarized and the other being transversely polarized, and _L_L for both 
J/'0s being transversely polarized. In these expressions, s, u are normal Mandelstam 
variables at the parton level. 
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Here, /j, and u are polarization indices of outgoing J/ ips, p and a are polarization indices of 
initial gluons, a and b stand for colors. 
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Here, rii — (1, — ps), = 0, and the polarized sum rule of s{XpJ is similar to that of £(Ap3). 
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